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Abstract. I apply the newly developed dynamical cluster approximation (DCA) 
to the calculation of the electron and phonon dispersions in the two dimensional 
Holstein model. In contrast to previous work, the DCA enables the effects 
of spatial fluctuations (non-local corrections) to be examined. Approximations 
neglecting and incorporating lowest-order vertex corrections are investigated. 
I calculate the phonon density of states, the renormalised phonon dispersion, 
the electron dispersion and electron spectral functions. I demonstrate how 
vertex corrections stabilise the solution, stopping a catastrophic softening of 
the (tt, tt) phonon mode. A kink in the electron dispersion is found in the 
normal state along the {{^, Q symmetry direction in both the vertex- and non- 
vertex-corrected theories for low phonon frequencies, corresponding directly to 
the renormalised phonon frequency at the (7r,0) point. This kink is accompanied 
by a sudden drop in the quasi-particle lifetime. Vertex and non-local corrections 
enhance the effects at large bare phonon frequencies. [PUBLISHED AS 
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1. Introduction 

Interest in the electron-phonon problem has recently been rekindled by the discovery 
of large electron-phonon couplings in the cuprate superconductors and the CMR 
manganites [2 [H Strong electron-phonon interactions outside the regime of 
traditional theories are also thought to be important in the bismuthates and in 
potassium doped BuckminsterfuUerene superconductors [3]. 

One of the best known approaches to the electron-phonon problem is that of 
Migdal and Eliashberg [H . Migdal showed that the effects of vertex corrections and 
momentum on the self-energy should be negligible if the phonon frequency is much 
less than the Fermi energy [5]. Migdal's method was developed further by Eliashberg 
to investigate the superconducting state [B . The physical content of Migdal's theorem 
is that there is a high probability that the most recently emitted phonon is the 
first phonon to be reabsorbed. This is true, provided excited states remain close 
to the Fermi-energy, i.e. phonon energies and couplings are small compared with the 
band width. Within Migdal theory, it is therefore acceptable to neglect processes 
in which the order of emission and reabsorption is changed. When these techniques 
were introduced, no materials had been found where the electron-phonon coupling 
was strong, and the phonon frequency large compared to the Fermi energy. For this 
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reason, the application of Migdal-Eliashberg (ME) theory was very successful, and 
remains highly regarded. Unfortunately, strong electron-phonon coupling and large 
phonon frequencies in the systems detailed in the first paragraph is likely to make 
them incompatible with the Migdal-Eliashberg approach, and further extensions to 
the theory seem necessary. The aim of this paper is to evaluate and discuss the effects 
of both vertex corrections and non-local fluctuations due to the DCA on the theory 
of coupled electron-phonon systems. 

Previous attempts to extend ME theory include the introduction of vertex 
corrections into the Eliashberg equations by Grabowski and Sham [7j, and an 
expansion to higher order in the Migdal parameter by Kostur and Mitrovic in order 
to investigate the 2D electron-phonon problem [5]. Grimaldi et al. generalised the 
Eliashberg equations to include momentum dependence and vertex corrections [S]. A 
discussion of the applicability of these and other approximations to the vertex function 
can be found in reference [10] . 

The importance of vertex corrections to electron-phonon problems within the local 
or dynamical mean- field approximation (DMFA), has been investigated by Miller et 
al. [H] and Deppeler et al. [12 . In these studies, the effects of both electrons 
on phonons and phonons on electrons are treated with equal importance, leading 
to a fully self-consistent theory. Miller et al. focus on situations where the effects 
of vertex corrections can be handled using a Coulomb pseudopotential extension to 
the basic ME theory, and where the explicit treatment of the vertex is essential. 
Deppeler et al. take the expansion to third order in the self-energy. The authors 
calculate physical observables, such as isotope coefficients, phonon spectral functions 
and superconducting transition temperatures. 

The current paper goes beyond the previous work by introducing a fully self- 
consistent momentum-dependent self-energy to the problem via the dynamical-cluster 
approximation (DCA). DCA is a new technique, which extends the DMFA by 
introducing short-range fluctuations in a controlled manner [13j . I apply the DCA 
to the calculation of electron and phonon dispersion curves in the Holstein model, 
which is one of the simplest non-trivial models of electron-phonon systems. Two 
approximations for the electron and phonon self energies are applied. The first 
neglects vertex corrections, but incorporates non-local fluctuations. The second 
incorporates lowest order vertex and non-local corrections. The vertex corrections 
allow the sequence of phonon absorption and emission to be reordered once, and 
therefore introduce exchange effects. At each stage, the DCA result is compared to 
the corresponding DMFA result. In principal, with parameters close to the Migdal 
regime, this enhancement of ME theory should be sufficient to correct the results, and 
the series allowing for additional reorderings should be convergent. However, such an 
approach may not work when the phonon frequency greatly exceeds the Fermi-energy, 
known as the instantaneous limit. 

This paper is organised as follows. In section[21 the DCA is introduced. In section 
[3]the Holstein model of electron-phonon interactions is described, and the perturbation 
theory and the full algorithm used in this work are detailed. In section U] the results 
are presented. They are divided into two parts. First the effects of non-local and 
vertex corrections on the phonon density of states and renormalised phonon frequency 
are discussed. Then the electron dispersion and spectral function are computed. A 
summary of the major findings of this research is given in section [S] 
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Figure 1. A schematic representation of the reciprocal-space coarse graining 

scheme for a 4 site DCA. Within the sliadcd areas, tlic sclf-cncigA' is assumed to 
be constant. There is a many to one mapping from the crosshatciicd areas to the 
points at the centre of tliose areas. Tlie coarse graining procedure corresponds to 
the mapping to a periodic cluster in real space, with spatial extent N^^^ ■ Also 
shown are the high symmetry points F, W and X, and lines connecting the high 
symmetry points. Since an infinite number of k states are involved in the many 
to one mapping, the approximation is in the thermodynamic limit. The DMFA 
corresponds to Nq = 1. 



n 




Figure 2. Diagrammatic representation of the current approximation. Series 

(a) represents the vertex-neglected theory which corresponds to the Migdal- 
Eliashberg approach. This is valid when there is a high probability that the last 
emitted phonon is the first to be reabsorbed, which is true if the phonon energy ojo 
and electron-phonon coupling U are small compared to the Fermi energy. Series 

(b) represents additional diagrams for the vertex corrected theory. The inclusion 
of the lowest order vertex correction allows the order of absorption and emission 
of phonons to be swapped once. In principal, for moderate phonon frequency and 
electron-phonon coupling, these additions to the theory should be sufficient to 
correct the results of the vertex-neglected theory. The phonon self energies are 
labelled with n, and T, denotes the electron self-energies. Lines represent the full 
electron Green's function and wavy lines the full phonon Green's function. 
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2. The dynamical cluster approximation 

The dynamical cluster approximation jl31 114j is a recent extension to the dynamical 
mean-field approximation. DMFA has been extensively applied to solving infinite- 
dimensional lattice models, where the DMFA formalism becomes exact [13 [TB]. It has 
also been used to approximate two and three dimensional systems, where it is known 
as the local approximation [T71[TT]. In DMFA, the self-energy is assumed to be local, 
and I](k, z) is replaced by its momentum-independent counterpart, S(2:). As such, 
the Green's function may be coarse grained across the Brillouin zone, 



where I?(e) is the non-interacting Fermion density of states (DOS). The chemical 
potential is represented by /i, and z may take the values iun or u + irj depending on 
whether one wants to work with the Matsubara or real axis Green's function (77 is a 
small value) . 

Application of the DMFA to one and two dimensional models is expected to 
give an incomplete description of the physics. The DMFA can be thought of as 
the first term of an expansion in l/unn, where n^^ is the coordination number |16j . 
For three dimensional systems, rinn ranges between 6 and 12, and it is often argued 
that additional terms give a minor correction. In contrasts additional terms in the 
l/rinn expansion are expected to be important in 1 and 2 dimensions. An example 
of significant differences between the one, two and three dimensional cases comes 
from isotropic-exchange spin-systems. In 3D, a finite transition temperature is found, 
consistent with the mean-field result (although the critical exponent is modified). 
However, significant non-local fluctuations in one and two dimensions are expected to 
reduce the Neel temperature to zero (Mermin- Wagner theorem), and the mean-field 
approach fails completely. 

The DCA follows a similar procedure to the DMFA. The Brillouin zone is divided 
up into Nc subzones consistent with the lattice symmetry (see figure [1]) . Within each 
of these zones, the self-energy is assumed to be momentum independent. This leads 
to the relation, 



where 2?i(e) is the non-interacting Fermion density of states for subzone i. The 
are the average k for each subzone, plotted as the large dots in figure [H This partial 
treatment of momentum dependence in the theory introduces non-local fluctuations 
with a characteristic length scale of Ny^. 

The combination of the self-energy and the coarse grained Green's function leads 
to a modified Dyson equation, 



The self-consistent condition is closed by calculating the self-energy from Qq{z,'K), 
which has the interpretation of the host Green's function of a cluster impurity model 
(sometimes known as the cluster excluded Green's function). Different approximations 
exist for the calculation of the self-energy, including quantum Monte-Carlo and 
perturbative approaches. 

Both dynamical-mean-field and dynamical-cluster approximations are calculated 
in the thermodynamic limit, as opposed to conventional finite sized techniques where 




(1) 




(2) 



Q^\^K^,z) - G-i(K„ z) - S(K„ z) 



(3) 
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the particle number is equivalent to the cluster size. In this work, the DCA is used, 
since it involves only a small amount of additional computing time per iteration, while 
rapid convergence in cluster size is expected. Estimates based on the convergence of 
the Hubbard problem (which has a relatively strong momentum dependence) suggest 
that the Hubbard model is well approximated by a cluster size of Nc = 64 [IS] . 

There are a number of methods that can be applied to speed up the DCA 
algorithm. Firstly, the symmetry of the problem may be taken into account through 
the appropriate planar point group pm3m [19]. In the large cluster limit, this means 
that only l/8th of the total momentum points need to be used in the calculation. 
Similar considerations can be used to dramatically reduce the number of terms in the 
momentum summations needed to calculate the self-energy terms. This is particularly 
important when vertex corrections are included, since the calculation scales as N^. 
The Matsubara frequency sums are, in principle, infinite. However, in this work they 
were truncated at a frequency where asymptotic behaviour is obeyed to a certain 
accuracy. All points above this frequency which were needed in the calculation of 
another quantity were computed from the asymptotic function, ensuring an accurate 
self-consistent procedure. The partial DOS are central to the DCA and were calculated 
using the analytic tetrahedron method to ensure very high accuracy f20^ . All code used 
in this study was written in object oriented C-I--I- to make the implementation as simple 
and reliable as possible. 

3. The Holstein model 

A simple, yet non-trivial, model of electron-phonon interactions treats phonons as 
individual nuclei vibrating in a static harmonic potential (representing the interaction 
between all nuclei) i.e. only one frequency ujq is considered. The phonons couple to 
the local electron density via a simple momentum-independent coupling constant g. 
The resulting Holstein Hamiltonian f2T| is written as, 

H = -J2 ^<y>-4cj- + J2 -l^) + J2 f^^^y^ + ^) (4) 

The first term in this Hamiltonian represents a tight binding model with hopping 
parameter t. Its Fourier transform takes the form gfe = — 2t ^^^j^ cos(fci). The second 
term connects the local ion displacement, Xi to the local electron density. Finally 
the last term can be identified as the bare phonon Hamiltonian. The creation and 
annihilation of electrons is represented by c|(ci), pi is the ion momentum and M the 
ion mass 

It is possible to find an expression for the effective interaction between electrons by 
integrating out phonon degrees of freedom [52]. In Matsubara space, this interaction 
is: 

m,..) = ^ (5) 

Here, u>s = ^-ksT are the Matsubara frequencies for bosons and s is an integer. A 
variable U = —g^jMuj^ is defined to represent the magnitude of the effective electron- 
electron coupling in the remainder of this paper. 

Interaction ([5]) has two important limits. When loq — > 0, the interaction is defined 
by a Kronecker b function. The DMFA of this limit was written down by Freericks 
et al, and solved by Millis et al. (H [23] . When phonon frequency and coupling are 
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small, Migdal's theorem applies. Migdal's approach allows vertex corrections to be 
neglected, and results in a momentum-independent theory corresponding to a local or 
dynamical mean-field solution. The theorem is exact when U — 0~ , ujo — ■ 

In the limit that lvq oo, the Holstein model maps onto an attractive Hubbard 
model. A particular feature of this mapping is that the second order direct and 
exchange diagrams of the expansion in the bare propagator develop the same functional 
form, and are of the same order of magnitude, while the first-order diagram becomes 
constant and can be absorbed into the chemical potential. Since the exchange 
diagram includes one vertex correction, it is clear that vertex corrections are extremely 
important in this limit. The vertex-corrected theory described in this paper has the 
appropriate weak coupling behaviour for large (jJq, unlike the vertex-neglected theory. 

In this paper, the perturbation theory shown in figure[2]is used, following directly 
from the free energy expansion of Baym and KadanofF [llj . The electron self-energy 
has two terms, I]me(w, K) neglects vertex corrections (figure [2J a)), and Svc('^,K) 
corresponds to the vertex corrected case (figure mb)). Hme('^,K) and Hvc('^,K) 
correspond to the equivalent phonon self energies. The diagrams translate as follows 
(note that in the following 4 equations, Q = (iw„, Q)) 



(K) = C/T^G(Q)i^(K-Q) (6) 



Q 



Hme(K) = -2C/r^G(Q)G(K + Q) (7) 



Q 



Evc(K) = {UTf G(Qi)G(Q2)G(K-Qi-Q2)i?(K-Qi)i?(Q2-Qi)(8) 

Ql:Q2 

Hvc(K) = -(C/r)2 J2 G(K)G(Qi)G(Q2)G(K-Qi+Q2)i?(Q2-Qi)(9) 
Qi,Q2 

The phonon propagator D{z, K) is calculated from, 

W^s,K) = , ^° — (10) 

i^i + ujo -H(iWs,K) 

and the Green's function from equation Q S = Eme + Eye- It is possible to 
extend equation [TO] to deal with a momentum dependent bare phonon dispersion, 
using a generalisation of the DMFA method in [SJ. When extending the formalism 
to investigate a momentum-dependent interaction term, the interaction term should 
be coarse grained as in [T^. All the ingredients now exist for the self-consistent loop. 
The algorithm is started with E,H = 0, then the DCA propagators, G and D are 
calculated from equations pU]) and H is then calculated followed by D and E. 
The loop is completed by calculating G. 

The double three-fold integration over momentum is the main barrier to 
performing vertex-corrected calculations, meaning that only relatively small clusters 
can be treated. Since the traditional phonon problem with uJo,U <C W (W is 
the bandwidth) has fluctuations which are almost momentum independent [5], it is 
expected that the DCA will have especially fast convergence in Nq for the parameter 
regime where ujo,U < W, and that calculations with relatively small cluster size 
accurately reflect the physics. In this respect, finite size calculations are likely to fail, 
and the application of DCA to this problem is essential. 
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Figure 3. The rcnormaliscd phonon DOS. The upper panel gives results for 
LOQ = 0.2, and the lower panel for LOa = 0.4. All four approximations arc shown 
on each panel. The vertical line represents the <5-function characteristic of the 
non-interacting DOS in both cases. The effect of electron phonon coupling is 
to shift spectral weight to lower frequencies, and to broaden the distribution, 
corresponding to a finite quasi-paxticle lifetime. Both panels show similar 
behaviour, but the effects are more pronounced for ujq = 0.4. The DMFA 
results show a simple distribution with a peak at the renormalised phonon 
frequency. Inclusion of vertex corrections in the DMFA leads to a reduction in the 
renormalisation, but the distribution remains similar. The inclusion of non-local 
corrections to the vertex-neglected theory causes an additional shift in the weight 
to very low frequencies. This shift is significant, because it causes a large increase 
in the renormalised electron-phonon coupling constant. The introduction of vertex 
corrections to the non-local theory causes the distribution to return to something 
similar to the original Migdal-Eliashberg result. This shows the importance of 
vertex corrections in stabilising the Migdal-Eliashberg theory. 



4. Results 



4-1- Phonon dispersion and DOS 

In this section, the phonon DOS and renormaUsed phonon frequency are calculated. A 
moderate couphng oiU = 0.3 was chosen, corresponding to approxiniatcly one third of 
the band width 1^ = 1 (i.e. U is slightly larger than t = 0.25). Two phonon frequencies 
of uuq — 0.2 and wq = 0.4 were chosen to demonstrate behaviour in the region where 
Migdal's theorem is (or ought to be) applicable, and behaviour in the region where 
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Figure 4. Renormalised phonon dispersion for the DCA. The solution from the 
DMFA results in a flat dispersion, and is therefore not shown here. The top row 
shows results for l^q = 0.2, and the bottom row for u)q = 0.4. Vertex neglected 
results are shown on the left hand side, and vertex corrected results on the right. 
As discussed in the text, the vertex corrections act to reduce the renormalisation 
of the phonon frequency. This can be seen to greatest effect in the = 0.2 
results, where the fully softened point in the vertex-neglected theory is no longer 
fully softened in the vertex corrected theory. The dispersion is relatively flat 
for the 0)0 = 0.4 case, which is reasonable since one expects a non-renormalised 
dispersion as one approaches the Hubbard limit {ujq ^ oo). 

vertex corrections are important. Computations were carried out with Nc — 1 and 
Nq = 64 in order to make comparisons between the DMFA and DCA result, with 
an additional computation for Nc = 16 in order to check convergence in cluster size. 
Although the results are not shown here, the Nc = 16 calculation was already close to 
convergence. A cutoff of 200 Matsubara points was used at a temperature of T = 0.02. 
The vertex neglected theory was calculated first, by iterating until convergence of 1 
part in 10^. The vertex corrections were then switched on, and iteration was continued 
until the required accuracy was reached. 

The total (momentum integrated) phonon density of states (cr(i-ij) — Im[D(cL;)]/7r) 
may be seen in figure [H Since it is difficult to choose a default model for the maximum 
entropy analytic continuation of Bosonic propagators [53], the analytic continuation 
in this subsection was carried out using the Fade approximant method of Vidberg and 
Serene [3^1. I also computed some results on the real axis for the vertex-neglected 
theory, and found no significant differences. 

The top panel in figure [3] shows results for luq = 0.2, and the bottom panel for 
Wo — 0-4- On each panel, results for all the approximations are shown to facilitate 
comparison (as is the case for most graphs in this section, exceptions will be noted with 
reasons). As is well documented for Migdal-Eliashberg theory, the effect of electron- 
phonon coupling is to shift spectral weight to lower frequencies, with the effect that the 
phonon frequency is renormalised to a smaller value (i.e. the phonon mode softens). 
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Figure 5. Electron dispersion of the Holstein model, shown for U = 0.3, = 0.2. 
Calculations for the DMFA and DCA with and without vertex corrections are 
shown. Note how the vertex corrections stabilise the kink. Some spurious 
dispersion is seen, and can be confirmed by looking at the spectral functions 
in figure (0 . However the weighting of these points is small, and they are marked 
with small dots for completeness. All four approximations show a kink at the 
renormalised phonon energy at the (vr, 0) point (the horizontal dotted line). The 
kink is especially well defined along the W-F line, where the 4 approaches agree 
to reasonable accuracy. It is possible that the kink is slightly below the phonon 
frequency for the DMFA calculations. A slight discrepancy is reasonable, since 
the renormalisation of the phonon frequency in DMFA is constant across the 
Brillouin zone. The kink is not well defined along the T-X-W direction, except 
for the DCA(VC) calculations, where a sudden discontinuity is seen. The point 
of this discontinuity does not coincide with a boundary to the coarse grained cell, 
so it is expected that it is not spurious. 



By way of comparison, the thick vertical hne represents the (5=function characteristic 
of the bare phonon DOS. 

The inclusion of vertex corrections into the DMFA calculation reduces the 
softening effect. This can be seen in both panels. The reduction in renormalisation 
occurs because the sign of live is opposite to that of IIme, i-e. the effect of particle 
exchange is to reduce correlation effects. 

The consequences of non-local corrections to the vertex-neglected theory are also 
shown in the figure. The effects are significant, since a proportion of the phonon 
spectral weight is shifted to very low energies. As the phonon frequency is increased, 
the relative size of the missing vertex corrections also increases. Therefore, the number 
of low frequency states rises. This very low energy behaviour is extremely important, 
since the small parameter for the Migdal-Eliashberg perturbation theory is normally 
written as A = "^9^ Jq <j{oj) duj/uj [27 . A is therefore significantly enhanced by the 
low energy modes, especially since the DOS at low energy has negative curvature. 
Examination of the partial phonon DOS for each coarse graining cell quickly shows 
that the main culprit is the (tt, tt) cell, which is fully softened. 

Vertex corrections to the non-local theory significantly reduce the renormalisation 
of the (tt, tt) point, causing the DOS to appear similar to that calculated from the local 
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Figure 6. As figure[5]but with U = 0.3, luq = 0.4. As before, vertex corrections 
stabilise the kink. Some spurious dispersion is seen, and can be confirmed by 
looking at the spectral functions in figure (0. The weighting of the spurious 
points is small, but they are marked with small points for completeness. Note 
that the kink energy is significantly lower than the appropriate renormalised 
phonon frequency at the (tt, 0) point for both of the DMFA calculations. The 
DCA(ME) dispersion bends smoothly (rather than showing a sudden kink), at 
an energy which is slightly lower than the phonon frequency. The DCA(VC) has 
kink behaviour along the W-F line in contrast to the other 3 approximations. 
Discontinuity in the dispersion along the F-X-W line can also be seen at this 
frequency. 



approximation. One expects this result, since the softening of the (tt, tt) mode increases 
the total magnitude of the vertex corrections within an iteration (by increasing the 
coupling constant A) reducing A over subsequent iterations. The phonon DOS finds 
an equilibrium between the softening of the vertex-neglected term and the effects of 
vertex corrections. This is consistent with experiments, where no full softening of the 
phonon frequency at the (tt, tt) point is observed for systems of optic phonons p, I28j. 
The effects of higher order vertex corrections are expected to be consistently smaller 
and smaller for this parameter range, since the lowest order vertex corrections were 
sufficient to stabilise A to a small value. Such an effect would presumably be irrelevant 
to systems of acoustic phonons, where the zone-centre mode is already zero. 

In order to elucidate the deficiencies in the vertex-neglected theory more 
completely, the renormalised phonon frequency in a coarse grained cell, f2(K), may 
be calculated by solving the equation, 

ujI^lu^ -LulRc[n{uj,K)] = (11) 

In figure m surface plots showing the renormalised phonon dispersion across the whole 
Brillouin zone are presented. Contours are plotted beneath the dispersion to highlight 
the resulting symmetry. Since the DMFA and bare phonon dispersion in the model are 
momentum independent, the renormalised phonon dispersion relating to the DMFA 
calculation of the Holstein model is constant across the Brillouin zone, and is therefore 
not shown. In a more realistic model, the dynamical effects of the coupling between 
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Figure 7. Spectral functions along the high symmetry directions for U = 0.3, 
LOQ = 0.2. Sharply peaked points close to the Fermi-energy (u) = 0) show well 
defined quasi-particles. These particles are suddenly damped at the renormalised 
phonon energy of the (7r,0) point (parallel dotted lines). However, at higher 
frequencies (F and W points), the damping is diminished in the vertex- neglected 
case. The vertex corrections predict that damping is maintained, even up to the 
largest energies. It is likely that the DMFA overestimates the phonon energy at 
(tt, 0). The full DC A with vertex corrections is necessary to obtain a very crisp 
change at the phonon frequency. 

nuclei would have to be taken into account, and one would start with a momentum 
dependent dispersion |29j . In order to apply DCA to such a problem, one would have 
to apply a coarse graining procedure to the phonon self-energy, in a similar spirit to 
|24| . In case a realistic model induces a momentum dependent coupling constant, then 
additional coarse graining could be carried out as in 14J. This would lead to a smooth 
dispersion as in the electron case. It is clear that the complete softening of the (tt, tt) 
mode that is seen in the vertex neglected theory (top left) is absent in the vertex 
corrected theory (top right). One of the main effects of non-local fluctuations can be 
seen here, since the DCA predicts that electron-phonon coupling induces a dispersion 
in the Holstein model (absent in the DMFA). As shown here, the dispersion induced 
from the coupling of electrons to a single Einstein mode is similar to the effects of 
inter-nucleus coupling in a non-interacting system. This means that it could be quite 
difficult to distinguish between the two effects in the analysis of experimental data. 

4-2. Electron dispersions and spectral functions 

The quasi-particle dispersion and lifetime are important physical quantities that can 
be measured by the ARPES technique. The renormalised dispersion can be calculated 
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Figure 8. As figure [71 but with U = 0.3, ojo = 0.4. In the case of the vertex 
corrected DCA, the quasi-particles are suddenly damped at the (tt, 0) phonon 
energy (shown by the parallel dotted lines). This sudden increase of damping 
is consistent with the experimental results in [2]. The vertex neglected theory 
consistently estimates the energy for the damping to begin at frequencies lower 
than the phonon frequency at the X point. As with figure Q, the vertex 
corrections maintain the damping up to the highest energies. 

by solving the equation, 

+ Ai - Ck - Re[E(w, K)] = (12) 

in analogy to the calculation of the phonon dispersion. Unlike the phonon calculations, 
the electron dispersion can be determined at all k points, owing to the many to one 
mapping involved in coarse graining the self-energy. Here K is the momentum of the 
coarse grained cell containing the momentum k. 

Figure [5] shows the electron dispersion of the Holstein model, for U — 0.3, 
uiQ = 0.2. The analytic continuation of the self energy was carried out using the 
maximum entropy method , since the normalisation of the spectral weight is known, 
and a reasonable default model can easily be constructed. The Fade approximant 
technique was not found to be very stable for the analytic continuation of self-energies. 
One reason for this may be that the self-energy has a high energy behaviour of 1/w, 
whereas the bosonic propagators have 1/w^ behaviour, and therefore the contribution 
to the approximant from higher frequencies is smaller. 

By examining the dispersions calculated for ujq = 0.2 it can be seen that the main 
feature of the electron dispersion is a "kink" at the energy corresponding to the phonon 
frequency at the X point (the phonon frequency is marked as a horizontal dotted line) . 
This is persistent for all the approximations. Vertex corrections stabilise the kink, and 
make it sharper. The kink is especially well defined along the W-T line, where the 
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Figure 9. Self energy of the Holstein model approximated using the DMFA, 
with and without vertex corrections. Since the self-energy of the DCA varies in a 
complicated manner throughout the Brillouin zone, it is not shown here. As the 
vertex corrections act against the Fock term, and have a much lower energy scale, 
the quasi-particle damping at the Fermi energy is suppressed. However, the zeroth 
order moment is uniquely defined by the Fock term, and so the total normalisation 
must be conserved. This means that weight is shifted toward higher energies, 
leading to the suppression of quasi-particle states at the top of the dispersion. 

4 approaches agree to reasonable accuracy. It is possible that the kink is located 
slightly below the renormalised phonon frequency in the dispersion calculated using 
the DMFA. Such a discrepancy is expected, since the renormalisation of the phonon 
frequency in DMFA is constant across the Brillouin zone, an assumption which is only 
true when phonon frequency and coupling are small as pointed out in section ITTl The 
kink is not well defined along the T-X-W direction, unless the dispersion is calculated 
using the DCA(VC) scheme, where a sudden discontinuity is seen. The location of 
this discontinuity does not coincide with the boundary to a coarse grained cell, and is 
therefore not a result of the DCA. I note that some spurious dispersion is seen (small 
points). However, the quasi-particle weight of these points is small, and they are only 
included for completeness. 

Figure [6] shows the electron dispersion for U = 0.3, = 0.4. Note that in the 
results from the DMFA calculations, the energy at which the kink occurs is significantly 
lower than the X-point phonon frequency. Results for the DCA(ME) bend smoothly, 
rather than kinking suddenly, at an energy scale, which may be slightly lower than 
the phonon frequency. Only the DCA(VC) has well defined kink behaviour along the 
W-T line. The discontinuity in the dispersion along the V-X-W line can also be seen 
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at this frequency. It is therefore quite clear that vertex corrections are significant for 
this set of parameters. 

In order to calculate spectral functions at general k, the coarse grained self energy 
is used in the following equation, 

k) = — — i — — (13) 

uj + irj + iJ. — 6]^ ~ 2j(tj, K) 

where K is the momentum point in the coarse graining cell which corresponds to k. 

In the calculation of the electron dispersion, some information is lost. For 
example, the quasi particle lifetime, which is related to the width of the peak in the 
spectral function cannot be seen. In order to demonstrate these effects, figure [7] shows 
a waterfall plot of spectral functions along the high symmetry directions for U = 0.3, 
loq ~ 0.2. The dotted lines running along the k axis represent the renormalised phonon 
frequency at the (tt, 0) point. 

Sharply peaked points can be seen close to the Fermi-energy (lu — 0), showing 
well defined quasi-particles. The quasi-particles are suddenly damped close to the 
renormalised phonon frequency corresponding to the (tt, 0) point, i.e. the peaks 
develop a width. At higher frequencies, the damping diminishes in the vertex neglected 
case and the peaks get higher and narrower. On the other hand, vertex corrections 
ensure that strong damping is maintained, even up to the largest energies. As 
discussed previously, it is thought that the DMFA overestimates the renormalised 
phonon frequency at (tt, 0). The full DCA with vertex corrections predicts a very crisp 
kink and reduction in the quasi-particle lifetime at the X-point phonon frequency. 

Figure [5] also shows spectral functions but now with the larger phonon frequency, 
luq = 0.4. As with the ujq — 0.2 results, the vertex-corrected DCA predicts a sudden 
damping of the quasi-particles at the (tt, 0) phonon energy. This sudden increase in 
damping connected with the kink is consistent with the photoemmision experiment 
in [2]. The vertex neglected theory predicts a much lower energy for the damping to 
begin. One can see that the kink predicted by the DMFA begins at a lower energy that 
the phonon frequency, and (for the vertex corrected DMFA) that the damping changes 
smoothly, rather than suddenly before the phonon energy, although full damping is 
not seen until above the phonon energy. 

The origin of the damping, and of the strong quasi-particles close to the Fermi- 
surface can be demonstrated by examining the self energy. Figure Q shows the 
self energy of the Holstein model approximated using the DMFA, with and without 
vertex corrections. The self-energy of the DCA is not shown, since it varies 
throughout the Brillouin zone. The form of the DCA self-energies are similar. As 
the vertex corrections act to reduce the self-energy at low energy scales, the self 
energy at the Fermi energy is suppressed, leading to well defined quasi-particles at 
low energies. However, the high frequency behaviour is defined by the first-order 
term in the self energy, and this can be related to the normalisation of the self-energy 
spectrum. Therefore the total normalisation is seen to be conserved, regardless of the 
approximation used. If the weight of the self-energy is reduced at low frequencies, this 
weight must be shifted toward higher energies, leading to the suppression of quasi- 
particle states at the top of the dispersion. 



5. Summary 



In this paper, I applied the newly developed dynamical cluster approximation to the 
calculation of the electron dispersion, induced phonon dispersion, spectral functions 
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and phonon DOS for the Holstein model. A theory neglecting vertex corrections 
(related to Migdal-Eliashberg theory) and another incorporating the lowest order 
vertex correction were investigated. In contrast to previous work, non-local corrections 
were included. The application of DCA to this problem was crucial, since the double 
three-fold (four-fold including Matsubara summations) integration involved with the 
vertex corrected sclf-encrgy of phonons and electrons means that only small clusters 
can be considered, precluding finite size calculations. 

The effect of non-local corrections to the vertex neglected theory of electron- 
phonon interactions is dramatic. Inclusion of non-local corrections to the theory 
results in a complete softening of the phonon modes at and close to the (tt, tt) point, 
which increases the effective coupling, A. Therefore, a first-principles approach to 
the electron-phonon problem in two dimensions which neglects the effects of vertex 
corrections, but keeps the effects of non-local corrections, can only be applied to a 
very small range of phonon frequencies and couplings. 

Vertex corrections act against the non- vertex corrected theory in both the phonon 
and electron self-energies. In the case of the phonon self-energy, which renormalises 
the phonon frequency, the correction acts to avert the complete softening of the mode. 
Therefore, A is stabilised. The theory with vertex corrections also acts to stabilise the 
kink in the electron dispersion at large phonon frequencies, and causes an enhancement 
in the kink. 

The results of this study have wide ranging implications for the application of 
Migdal-Eliashberg theory to the electron-phonon problem. Dispersion can be induced 
in the phonon modes, making it difficult to estimate the true coupling from fits 
to experimental data. Vertex corrections act to return the system to the Migdal 
parameter regime. Therefore, it is not surprising that the Migdal-Eliashberg theory 
has been so successful in describing the properties of a wide range of superconductors. 

There is a large scope for extension to this work. The most direct and 
straightforward extensions would be to investigate the superconducting phase, the 
phase diagram and the 3D problem. It would also be interesting to study the effects of 
even higher-order vertex corrections. This is not practical within perturbation theory, 
since an implementation would involve would involve many 3-fold integrations. An 
alternative method for including vertex corrections would be to employ a quantum 
Monte-Carlo simulation. Such work is currently in progress. 
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